Exact c-number Representation of Non-Markovian Quantum Dissipation 
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The reduced dynamics of a quantum system interacting witli a linear heat bath finds an exact representation 
in terms of a stochastic Schrodinger equation. All memory effects of the reservoir are transformed into noise 
correlations and mean-field friction. The classical limit of the resulting stochastic dynamics is shown to be 
a generalized Langevin equation, and conventional quantum state diffusion is recovered in the Born-Markov 
approximation. The non-Markovian exact dynamics, valid at arbitrary temperature and damping strength, is 
exemplified by an application to the dissipative two-state system. 
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Irreversible quantum processes are important in almost ev- 
ery field of condensed-matter physics and chemistry. In phe- 
nomena as different as diffusion of light particles in solids 
and light harvesting in biological systems, the dissipation of 
energy, the destruction of phase coherence, and the genera- 
tion of entropy play a key role. The question how such pro- 
cesses can be described by the stochastic propagation of quan- 
tum states has recently attracted increasing attention 
For classical systems with linear dissipation, Langevin equa- 
tions provide a theoretical (and numerical) tool to accurately 
describe the interaction of a system with a complex thermal 
reservoir in comparably simple terms of stochastic forces and 
memory friction. For open quantum systems, no such sim- 
ple and exact approach in terms of an equation of motion has 
been established so far The dynamical quantity of interest is 
the reduced density matrix, obtained by tracing out the reser- 
voir degrees of freedom from the dynamics. Significant in- 
formation about correlations between system and reservoir is 
lost in this operation, hence no deterministic differential equa- 
tion of motion can be derived for the reduced density matrix 
without approximations. Traditional approaches treat these 
correlations perturbatively, yielding approximate equations of 
motion such as Redfield or Master equations. While these 
have been used very successfully in the fields of quantum 
optics and magnetic resonance, there are many condensed- 
matter problems where they are qualitatively wrong. Large 
coupling constants and long correlation timescales both need 
to be treated non-perturbatively. This can be accomplished in 
a formally exact manner by path integrals for open quantum 
systems [^-^, but an analytic evaluation of these functional 
integrals is usually restricted to special cases and approxima- 
tions, e.g., the semiclassical limit. Furthermore, present exact 
numerical techniques such as Quantum Monte Carlo (QMC) 
methods | ]lC| , pl] ] have to cope with a sign problem arising 
from the fact that probability amplitudes with varying phases 
rather than probabilities have to be added in a quantum com- 
putation. 



Quantum State diffusion (QSD) has been established 
as an alternative theory of quantum dissipation in the pertur- 
bative regime. Through the stochastic propagation of pure 
quantum states it yields an intuitively appealing picture of the 
behavior of individual quantum trajectories in open systems 
and permits effective numerical calculations. Generalizing 
this formalism to the case of non-perturbative, non-Markovian 
dynamics holds the promise to overcome most of the above- 
mentioned limitations of currently known techniques. This 
approach has recently been taken by Diosi, Strunz, and Gisin 
[Q], but at the price of incorporating non-Markovian retarda- 
tion effects in a memory functional. In practice this means that 
a general solution of the resulting equation of motion, even 
numerically, is almost as elusive as that of the underlying path 
integral. 

In this Letter we derive and discuss stochastic Schrodinger 
equations for open systems which allow the treatment of linear 
dissipation of arbitrary strength and correlation time scales. 
Without approximation, all effects of the system-reservoir in- 
teraction are re-cast in the form of c-number noise and friction 
forces with a suggestive physical interpretation: The ensem- 
ble of phase-space points described by a classical Langevin 
equation is generalized to a stochastic ensemble of quantum 
states. Both descriptions are linked by a correspondence prin- 
ciple for open systems. 

Any general theory of quantum dissipation has to start from 
an open system embedded in a large reservoir system, whose 
degrees of freedom are treated fully quantum mechanically, 
but can later be eliminated from the dynamics. The Hamilto- 
nian of such a model consists of system and reservoir terms 
and an interaction potential. 



H = i7o + Hji + Hi . 



(1) 



Let us consider a thermal coiTelation function 

{AB{t)) = trjc-'^^^e^^'Be-^^*} /tie-^" , (2) 
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where A and B are system operators. The time evolution 
operators can be joined with the "imaginary-time" propaga- 
tor e~^^ to form a propagator with complex time argument 
which describes the time evolution of a quantum system along 
the contour C depicted in Fig. |l|. This time evolution can be 
expressed by a path integral over the configuration space of 
the system 

{ABit)) oc J P[q]e^^«['ilj[q](qo+|A|qo-)(qt+|i?|qt-), 

(3) 

where a normalization factor has been suppressed for simplic- 
ity. The action functional 



^o[q] =^clTi(q,q) 



(4) 



depends in the usual way on the classical Lagrange function. 
The measure dr is the complex time differential on C, and 
q = . In the generalized influence functional 

I[ci] = (tcCxp(^^ J^dT^Hj{q)^^ (5) 

the thermal average of the contour-ordered exponential is 
taken with respect to the unperturbed reservoir This exact 
result can be simplified [p"2p|] by performing a cumulant ex- 
pansion of the expectation value in I[q\. The most commonly 
discussed case is that of linear dissipation, equivalent to trun- 
cating the cumulant expansion after the second order This 
procedure is exact for a model reservoir of harmonic oscilla- 
tors, but also applies to any microscopic model whose collec- 
tive response to a perturbation Hi is dominated by its linear 
term. 

In the following we restrict ourselves to linear dissipation, 
for which /[q] is a Gaussian functional determined by the cor- 
relation function {Hj{q{T))Hi{q{T'))). After an expansion 
of Hi{q) in terms of an arbitrary set of reservoir operators 



Hi = Y,fM)x, 



(6) 



the correlation matrix Ljk{T — t') = {Xj{T)Xk{T')) p, char- 
acterizing the isolated reservoir, completely determines the in- 
fluence functional 

/[q]=exp( / dr / dr'/,(q(r)) 



L,k{r^r')M<l{T')) 



(7) 



Here denotes the order relation induced by the contour 
orientation. 

The aim of the present work is to transform the exact de- 
scription of open quantum systems by the influence func- 
tional technique into an equivalent stochastic propagation of 



the quantum states of the open system. We first show that the 
influence functional /[q] can be constructed as a noise average 
of the form 



/[q] = (exp(^^ / drz,(r)/,(q(r)) 



(8) 



w 



where the average is taken with a Gaussian probability mea- 
sure If the stochastic co variance matrix of the 
complex- valued noise coordinates zj (r) matches the quantum 
correlation matrix of the coupling operators Xj, 



\Ljk{T - r'), T ^ r' 

(Zj(T)zfc(T ))w = { 

\Ljk{T' -t), t r', 



(9) 



Eqs. and ^ constitute a formal identity between Gaus- 
sian functional. Noting that Eq. does not fully deter- 
mine the noise statistics, we find that correlations of the type 
{zj{T)z1.{T')) can be chosen such that the covariance ma- 
trix of all real-valued noise components represents a posi- 
tive quadratic form, i.e., Gaussian noise {zj} with the desired 
properties exists. 

For each sample of the noise, the propagation is now gov- 
erned by a time-local action functional, i.e., the path-integral 
dynamics can be translated into an equivalent Schrodinger 
equation. The propagation along the contour C is governed 
by the Hamiltonian H^ and a stochastic potential term. 



\h\i>) = Ho\i>) -y,z,{t) fMm- 



(10) 



This remarkably simple equation shows that linear dissipation 
can be described exactly by a linear QSD theory containing 
no memory terms. 

In order to compare this new finding with previous re- 
sults, we need to choose the Feynman- Vernon influence func- 
tional [H], for which the noise vanishes on the imaginary 
axis. For simplicity, we discuss a one-dimensional coupling 
Hi = —qx, where x is a reservoir operator with correlation 
function L{t — r'), associated with a single noise coordinate 
z{t). The propagation along the two real-time segments of the 
contour can be re-stated in the form of a stochastic Liouville 
equation for the reduced density matrix. 



= [^^0, p] - ziqp + z^pq, 



(11) 



where zi{t) — z{t) and Z2{t) = z{t — ifif3), and where 
the sample p of the reduced density matrix is separable, p — 
\ipi){ip2\- Hence Eq. (11 ) is just a compact notation for two 
stochastic Schrodinger equations for and \ip2)- 

The noise forces may be represented as the sum of statisti- 
cally independent terms, zi.2 = z + vi,2- The terms vi and V2 
have identical statistics, but are uncorrelated. Performing an 
average over wi 2 takes us back to the partial decomposition 
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of the influence functional given by Diosi, Strunz and Gisin 
10]. Using the same representation, the stochastic Schrodinger 
equation of Markovian QSD can be recovered. The same Um- 
its and approximations used to derive Lindblad-type Master 
equations and their QSD counterparts from a system-reservoir 
model, i.e., weak damping and time coarse-graining, can be 
applied to Eq. (^Tj). Coarse-graining is done at an intermedi- 
ate timescale which is short compared to the relaxation time 
of the damped system but long compared to its oscillation pe- 
riods and the relaxation times of the reservoir. In this approxi- 
mation the dynamics is only sensitive to segments of the noise 
spectrum of z which are centered around the natural transi- 
tion frequencies of the undamped system. In the Born-Markov 
limit the coarse-graining time scales to zero, and each of these 
segments naturally transforms into a distinct Markovian noise 
force coupled to a Lindblad operator associated with the re- 
spective quantum transition. Finally, the work of Stockburger 
and Mak [|]] using real noise and earlier work concentrat- 
ing on strictly ohmic damping [ pj| ] can be recovered from 
the present decomposition of the Feynman- Vernon influence 
functional by applying the stochastic decomposition only to 
the real part of L(t — t') while explicitly evaluating memory 
effects induced by the imaginary part. 

As in the linear versions of conventional QSD, tr p is not 
conserved by the time evolution of Eqs. ([l^) and (pi]). A 
transparent physical interpretation of the noise forces and 
stochastic samples arises in a modified nonlinear theory for 
which the trace of the reduced density matrix is preserved not 
only in the stochastic average, but for each sample. 

It will be advantageous to introduce new noise variables 
^ = (zi + 22)72 and fiv — z\ — Z2, with resulting covariances 

{m^it'))w ^ReLit-t') (12) 
{^{t)iy{t'))w = {2i/n)e{t - t') lmL{t - t') 

= -ixB{t - t') (13) 

{u{t)u{t'))w = 0, (14) 

where XR is the response function of the reservoir. The noise 
force ^(t) can be chosen real, and this case will be discussed 
here for simplicity. For comparison with phenomenological 
friction models, one conventionally augments the coupling 
Hamiltonian by a counter-term in order to make the coupling 
translationally invariant. The additional term is a quadratic 
potential modification iiq^ /2 with fi = J^^ dtxR{t), which 
eliminates the static response of the reservoir. Eq. ( |ll| ) is thus 
transformed into 



ihp = [Ho,p] + -C[q,p] - ^''{q,p}- 



(15) 



For the normalized density matrix sample p = p/ tr p, this 
translates into the quasiUnear equation of motion 

in'p^[Ho,p] + ^[q^,p]~C[q,p]-^l^{q~q,p} (16) 

with ^ = tr qp. When averaging the normalized density 
matrix p, the factor 



trp = cxp ^i J dt'q(t')iy{t')^ 



(17) 



needs to be incorporated into the integration measure. Be- 
cause the original integration measure j/, z^*] is Gaus- 
sian, the new measure Wt[£,,i^,h'*] including the factor tr p 
can be rewritten by "completing the square" in the exponent 
of the Gaussian functional, making the exponent a quadratic 
form of shifted variables, and yielding VFt[^, i^, i^*] = 

with 



(18) 
(19) 



^t{t') = m+ I dsxRit' ~s)q{s) 



ly^it') = I,* it') -i ds {u*{t')u{s))wq{s) (20) 
Jo 

Using W['Ct, v'^] as integration measure allows us to in- 
terpret the shifted variables as noise to be used for stochastic 
integration. The Jacobian determinant of the variable change 
is unity: Because vt = v, and because q is independent of v*, 
it simplifies to 
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(21) 



Using the fact that both the bath response and v\ are 
causal, one finds that |(5^t/(5^| is represented by a triangular 
matrix with unit diagonal elements, i. e., |(5^t/(5^| = 1. 

After the variable change (with subscripts t now dropped) 
the equation of motion is rewritten as 

1^/5= [-ffo,/5] -^[(7,/5] +m-f{t)q{Qi)[q,p] 

+ mf dt'j{t^t')q{t')[q,p] 
Jo 



+ 2^i<l-Q) ,P] - 2^{'?-9,p}• 



(22) 



Here we have integrated by parts, using the relation mj{t) = 
— X7? (i) between the reservoir response function and the fric- 
tion kernel 7(i). Remarkably, the mere choice of a con- 
straint for tr p has taken us from the abstract mathematical 
decomposition (||) to an equation of motion with a number of 
terms looking quite familiar — as in classical dissipation, the 
dynamics is governed by c-number force terms representing 
noise and a mean-field version of classical friction. Additional 
terms ascertain that the noise-averaged p contains the effect of 
all system-reservoir correlations. 

Further examining this analogy, we find a "sample-by- 
sample" correspondence between quantum mechanical and 
classical stochastic dynamics. The time evolution of Eq. ( p2[ ) 
can be described in a reduced Heisenberg picture, which is de- 
fined by interpreting expectation values of the form tr{Atp} 
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as a scalar product {A, p) of two operators. This allows the 
introduction of Heisenberg operators A{t) through the propa- 
gating superoperator and its adjoint, 

{A,p{t)) = [AM{t)M = iUHt)A,po) = iAit),po). (23) 



In the classical limit, terms proportional to q—qin Eq. ( |22| ) can 
be neglected, and the equations of motion governing the evo- 
lution of Heisenberg operators A{t) become those of a classi- 
cal stochastic system. For a potential model, these reduce [ p^ 
to the generalized Langevin equation 

mm = -V'iq{t)) -m f dt'jit - t')q{t') + m ■ (24) 

The limit fi ^ Q turns our stochastic ensemble of quantum 
states gradually into an ensemble of phase-space points with 
classical thermal noise and memory friction, i.e., the quantum 
dynamics defined in Eq. ( |2^ obeys a stochastic correspon- 
dence principle. 

Fig. 2 shows numerical results for the symmetric spin- 
boson system [p3[], where q is the Pauli matrix and 
Ho = (?iA/2)(Tj;, compared to a QMC algorithm which di- 
rectly samples a path integral [p^. For the coherent dynam- 
ics shown in Fig. 2a, the QMC method does not fully con- 
verge even if run longer than our simulation. Due to an only 
marginally ergodic Metropolis random walk, the accuracy of 
the QMC method degrades significantly with increasing t. 
Problems of this kind never appear in our algorithm because 
our samples are statistically independent by construction. The 
incoherent dynamics shown in Fig. 2b, a forte of the QMC 
method, is treated equally well by both approaches. 

We have derived exact stochastic equations of motion for 
the reduced density matrix in linear form ( [lo| ) and trace- 
conserving form (P2|). In the latter, a close correspondence 
with classical Langevin equations is evident. Examples of nu- 
merical results show that the stochastic equations of motion 
presented here provide a promising alternative to currently 
available exact methods based directly on the path integral 
representation. Given the freedom of choice for the noise co- 
variances left undetermined by the physics of quantum dissi- 
pation, future improved numerical methods are to be antici- 
pated. 
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FIG. 1. Contour for propagators in thermal correlation functions. 
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FIG. 2. Stochastic and quantum Monte Carlo simulations of a 
symmetric spin-boson system with XrO-) c< texp{—ujct), dissipa- 
tion constant [^] Q = 0.1. (a) Coherent dynamics at zero temper- 
ature. (Five repeated simulations of Eq. (llSp are shown to indicate 
statistical errors), (b) Incoherent dynamics at T = 5hA/kB, simu- 
lated using Eq. (|22|). 
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